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Guided atom laser: transverse mode quality and longitudinal momentum distribution 
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We analyze the outcoupling of a matter wave into a guide by a time-dependent spilling of the atoms 
from an initially trapped Bose-Einstein condensate. This process yields intrinsically a breakdown 
of the adiabatic condition that triggers the outcoupling of the wave function. Our analysis of the 
time-dependent engineering and manipulation of condensates in momentum space in this context 
enables to work out the limits due to interactions in the mode quality of a guided atom laser. This 
study is consistent with recent experimental observations of low transverse excitations of guided 
atom lasers and suggests (i) an optimal strategy to realize such quasi-monomode guided atom lasers 
with, in addition, the lowest possible longitudinal velocity dispersion, or alternatively (ii) a strategy 
for engineering the atomic flux of the atom laser. 

PACS numbers: 03.75.Pp,37.25.-|-k, 41.85.Ew, 42.60.Jf 
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Well-controlled coherent atom sources such as atom 
lasers are desirable for atom interferometry [1[ and can 
be used as quantum probes of complex externalpoten- 
tials to study quantum transport phenomena [2|, |3| . The 
first generation of atom lasers were obtained by outcou- 
pling atoms from a trapped Bose-Einstein condensate 
(BEG) to free space [J|. Different aspects of these free 
falling atom laser have been extensively studied: the op- 
timization and the characteristic of the output flux (5|-l8j , 
the longitudinal velocity dispersion Q and the transverse 
mode quality [13, Ell- 

However, the de Broglie wavelength of a free falling 
atom lasers decreases rapidly because of the acceleration 
under gravity. To circumvent this problem a second gen- 
eration of atom lasers involves atom outcoupling in an 
horizontal optical guide |12h15| . In this way, a mean de 
Broglie wavelength on the order of one micrometer can 
be achieved over macroscopic distances of a few mm. 

Similarly to free falling atom lasers, different outcou- 
pling schemes have been demonstrated. In 12], the trap 
is magnetic and atoms are outcoupled by radio- frequency 
spin flip. In [Ij], atoms are extracted from a BEG held 
in a crossed dipole trap by applying an increasing mag- 
netic field gradient along the guide direction. Alterna- 
tively, optical outcoupling was demonstrated and consists 
in lowering progressivelly the intensity of one of the beam 
of the crossed dipole trapped to favor the departure of 
atoms in the other beam of the trap [1Jt3 ■ 

The reproducible preparation of a well-controlled prop- 
agating matter-wave is always a delicate task and is gov- 
erned by the control of the extraction process. The mean 
number of modes occupied in the transverse direction 
can be easily infered from time-of-flight expansion for 
alkali atoms jl3 - [l4| or by a detector based on a mi- 
crochannel plate (MGP) for metastable helium atoms 
[16| . Gertain schemes prove the possibility to realize a 
quasi-monomode guided atom laser (GAL) while other 
schemes remain multimode with a small number of pop- 
ulated transverse modes [ij, H, Il3 ■ Technical noise can 
affect dramatically the quality of the transverse mode of 
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FIG. 1: Sketch of the optical outcoupling procedure. The 
intensity of the beam that ensures the longitudinal trapping 
is reduced progressively as a function of time which favors the 
spilling of the wave function outside the trap region. 



GAL. For instance, magnetic noise may lead to a limit in 
the control of rf outcoupling schemes jl7| . However, even 
if one overcomes such experimental difficulties, one may 
wonder, for a given outcoupling process, to which extent 
there exists a fundamental limit introduced by the out- 
couling mechanisms. A related question is to know what 
is the best strategy to realize an optimal outoucoupling 
for both longitudinal and transverse degrees of freedom. 
Indeed the longitudinal dispersion velocity governs the 
longitudinal coherence length of the atom laser. In all 
GAL experiments, there is a lack of information about 
this quantity and there is a need for methods allowing 
one to measure and engineer it |18| . Indeed narrow mo- 
mentum distribution are required to obtain good fringe 
contrast in atom interferometry experiments |l|. 

In this article, we focus on the study of the optical 
outcoupling scheme similar to the one implemented in 
14 - 16] . Its principle is sketched in Fig.[TJ Atoms are out- 
coupled from the trap by lowering its depth as a function 
of time. When the trap depth decreases, the strength 



of the trap also decreases and the minimum of the trap- 
ping potential is slightly displaced. At some point, even 
if the lowering of the trap depth is very slow, the adia- 
batic evolution breaks down since the timescale associ- 
ated with the instantaneous trap diverges when the trap 
depth is reduced. This breakdown triggers the outcou- 
pling of the wave function. Strictly speaking, this phe- 
nomenon rules out a priori any attempt based on a naive 
entropic point of view according to which the ground 
state wavefunction of the initially trapped atoms can be 
mapped on the ground state of the guide. To analyze this 
time-dependent problem, we mainly resort on numerical 
simulations. 

The paper is organized as follows. In Sec. I, we discuss 
the impact of the breakdown of adiabaticity on the lon- 
gitudinal momentum dispersion of the atomic wavefunc- 
tion in a one-dimensional geometry and in the absence 
of interactions. In Sec. II we take also into account the 
interactions within the mean field theory. In Sec. Ill, we 
discuss the possibility of engineering the atomic flux us- 
ing a simple model. In Sec. IV, we use a two-dimensional 
geometry to account for the transverse modes in the guide 
and include the interactions. Our results explain the 
quality of the transverse mode obtained experimentally 
[l3 | and provide a new strategy to reduce at best the 
longitudinal velocity dispersion requirements while pre- 
serving a very good transverse mode quality. 



I. ID OUTCOUPLING OF AN IDEAL BOSE 
GAS 

The potential experienced by the atoms is chosen in 
the form 



U{x, t) = ^ra-ix - [/o/(t)e-2^'/" 



(1) 



The function /(i) characterizes the time dependence of 
the trap depth and tends to zero when time tends to in- 
finity to ensure the outcoupling. The other parameters 
are chosen to reproduce the typical experimental values: 
the mass m is the one of rubidium-87 atom, the slope 
(first term of the r.h.s. of Eq. ([1])) is charaterized by an 
acceleration 7, and the angular frequency at the bottom 
of the trap is wq — 2\/TJZ/'Wo = 27r x 150 s~^ (depth 
Uo/kB = 3.7 fiK and waist wq — 40 /im). The param- 
eters are such that the time associated with the tunnel 
effect remains always extremely large compared to all 
other timescales of the problem. The initial state of the 
wave function is here the ground state wave function as- 
sociated with the harmonic oscillator potential obtained 
by the expansion of U{x, t) to the second order in x about 
the minimum x = xq, dxU{x,0)\xo — (see Fig.[T]). 

We evolve the condensate wave function using a split- 
operator fast Fourier transform method [I3. The evolu- 
tion of the velocity dispersion in the course of the out- 
coupling processes is plotted in Fig. [2] for an exponential 
time function f(t) = e"*/"^ and different characteristic 
times r. Such an exponential form has been used in 
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FIG. 2: Ratio between the longitudinal velocity dispersion of 
an ideal Bose-Einstein wavefunction as a function of time in 
the course of the outcoupling process and the initial veloc- 
ity dispersion Ai;(0) — 0.58 mni.s"^ for different outcoupling 
characteristic times t — ms, t — 1 ms, 5 ms and 100 ms. 
The time to for which the potential does not have anymore a 
minimum is io — 4r. 



[13| . An abrupt switch off of the trap (r — 0) enables 
to keep the velocity dispersion constant. Indeed, in this 
case, the whole wave function experiences only the lin- 
ear potential slope and it is a well known result, used in 
time-of-flight experiments, that the velocity dispersion of 
a wave function is not affected when atoms are subjected 
to a constant acceleration field. When the characteris- 
tic time, T, is non-zero, the dispersion velocity starts to 
decrease. This is due to the beneficial effect of a quasi- 
adiabatic opening of the trap which increases the size 
of the cloud and decreases correspondingly the momen- 
tum dispersion since the ground state is a minimal state 
{AxqApo = ^/2). The maximum gain in the reduction 
of the momentum dispersion is obtained sd t ^ to for 
T ^ 5 ms where to is defined by 
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for the chosen exponential form of the time function f(t) 
and corresponds to the time at which the potential does 
not have anymore a local minimum. 

For longer time as used in experiments |14l . |l5| (see 
Fig. [2] for r = 100 ms), three different phases can be 
distinguished. First, the velocity dispersion decreases as 
a result of the adiabatic opening of the trap. Second, a 
rapid increase occurs. This corresponds to the outcou- 
pling stage. The large velocity dispersion originates from 
the fact that a part of the wave function is accelerated 
on the slope while the other part is still in the trap. This 
results in a stretching of the wave function in momentum 
space. The last stage corresponds to the propagation of 
the packet on the slope (constant momentum dispersion). 
The minimum value of the achievable velocity dispersion 
depends on the slope 7. In the limit of a vanishing slope, 
the lower bound is given by the Heisenberg limit h/2mwQ. 

We conclude that, in the absence of interactions, an 
optimal strategy to minimize the longitudinal velocity 
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FIG. 3: Wigner function of the wave function at different 
times t/r = 0, 4.5 and 7 for f{t) = e"*/^ with r = 10 ms. The 
wave function is evolved according to the ID Gross Pitaevskii 
wave function with an initial condition corresponding to the 
Thomas Fermi regime {Ngi = 5.ih^'^ujQ m,~^'^). 



dispersion, and referred in the following as the two stage 
strategy, consists in decompressing progressively the trap 
up to the point where outcoupling occurs (i.e. for i ~ io) 
and then to switch it off abruptly so to freeze the veloc- 
ity dispersion. In this way, it is possible to decrease the 
longitudinal velocity of the trapped atoms by a factor 
slightly larger than two for 7 = 0.2 m.s^^. This strat- 
egy does not depend on the specific shape of the function 
fit). It is worth mentionning that the adiabatic expan- 
sion stage can be speed up using shortcut methods as 
recently demonstrated |20l - l23l |. 



II. ID OUTCOUPLING OF AN INTERACTING 
BOSE-EINSTEIN CONDENSATE 



The description of a condensate in the presence of re- 
pulsive interactions is made through an extra mean field 
term in the Schrodinger equation. The corresponding ID 
Gross Pitaevskii (GP) reads 



ih 



~dt 



2m dx^ 



+ C/(a;,i)^ + gi|V|V, (3) 



where the wave function ^ is normalized to the total 
number of particles. At i = 0, the BEG wave function is 
taken to be a stationary solution of the time independent 
GP equation. It is found by evolving the time-dependent 
GP equation in imaginary time i — > ii in the potential 
dD) with 7 = (no slope) but taking into account the 
displacement of the minimum xi^ induced by the slope. 
This trick facilitates the numerical convergence towards 
the stationary state. We have chosen the strength pa- 
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FIG. 4: Ratio between the longitudinal velocity dispersion of 
a ID wave packet and the initial velocity dispersion At;(0) = 
0.25 mm.s~^ as a function of time in the course of the out- 
coupling process for (i) a sudden switch off of the trap (rep- 
resented in units r = 1 ms), and (ii) for different outcoupling 
characteristic time: r =1 ms, r =5 ms, r =100 ms. Gray 
solid lines: numerical simulation. Dashed lines: results from 
the analytical model in momentum space for the opening of 
the trap detailed in Appendix A. The breakdown of the good 
agreement between simulation and the model occurs when the 
spilling of atoms starts. The wave packet describes the wave 
function of a condensate initially in the Thomas Fermi regime 
with an interaction parameter Ngi = 53h^'^u!Q m~^'^. 



density profile is close to an inverted parabola: 



ni{x,t = 0) = 



37V 



[R'l? 



(4) 



condensate is well in the Thomas Fermi regime and the 



where the Thomas Fermi radius reads R\ = 
{iN gi/2mujQY/^ ~ 3.8 /im with our parameters [24| . 
With these parameters close to those used in experiments 
[lJill3; ^^(0) = 0.25 mm.s"^. The interaction energy 
per particle reads E'^\R\) = 2^;/5 = 37Vgi/10i?? where 
/i]* is the chemical potential. 

Figure [3] shows the Wigner function in a phase space 
representation of the wave function at three different 
times: in the trap before the outcoupling, just after the 
outcoupling procedure, and in the course of the expan- 
sion. A space- velocity correlation develops as a function 
of time. As a result, the analysis in position space of 
the interaction of an atom laser with an extra potential 
provides directly the behavior for the different class of 
velocity present in the wave packet. This feature is in- 
teresting in view of using atom lasers as quantum probes 

The evolution of the momentum dispersion of the wave 
function as a function of time is plotted in Fig. H] for dif- 
ferent characteristic times r. The momentum dispersion 
increases at short time t < t (before any outcoupling) for 
low values of the characteristic time r = 0, 1 ms but de- 
creases for longer characteristic times r = 100 ms. This 



means that two contributions with opposite effect are in- 
volved in the dynamics of the momentum dispersion. At 
an intermediate characteristic time r = 5 nis, this trade- 
off yields a transient oscillatory behavior (see Fig. U). 
The increase is well captured by the standard analysis 
based on scaling parameters of a Thomas Fermi cloud. 
It is due to the conversion of interaction energy into ki- 
netic energy. The decrease is attributed to the adiabatic 
expansion and corresponds to the decrease of the proper 
kinetic energy. The dynamics of the trade-off in between 
those two contributions is discussed in Appendix A using 
an extension of the standard scaling ansatz method. The 
corresponding predictions based on an harmonic approx- 
imation of the bottom of the trap are shown as dashed 
lines in Fig. |4l They capture very well the evolution of 
the momentum dispersion in all range of values for the 
characteristic time r, and are valid during the opening of 
the trap and before the outcoupling of atoms. 

For short characteristic times (r <1 ms), we obtain 
a poor gain (on the order of 10 %) in the asymptotic 
longitudinal dispersion velocity compared to its value for 
an abrupt switch off of the trap. For r ~ 5 ms, a gain 
up to 50 % can be obtained as a result of the beneficial 
effect of the adiabatic opening of the trap. As previously, 
the two stage strategy is still efficient to get the smallest 
possible longitudinal velocity dispersion. However, such 
a strategy prohibits to engineer the atomic fiux (see next 
section). 

When the strength of the repulsive interactions is 
increased, the initial size of the condensate increases 

1/3 

{R oc g^ in the Thomas Fermi regime) and the mo- 
mentum dispersion decreases compared to its value in 
the absence of interactions. However and in contrast 
with the ideal Bose gas case, the presence of repulsive 
interactions yields to an increase of the asymptotic ve- 
locity dispersion Au(oo) compared to its initial value 
Ai;(0) in all schemes, i.e. whatever is the characteris- 
tic time r. This is due to both the transfer of the in- 
teraction energy into kinetic energy in the course of the 
expansion of the wave packet (dominant process at low 
r values) and the stretching effect of the wavefunction in 
the course of its outcoupling (dominant process at large 
r values). For instance, assuming t — 5 ms, we find 
Av(t — > oo;gi = gf^) ~ 5.5Av{t -^ oo; gi — 0). Interest- 
ingly, the timescale to reach the asymptotic value turns 
out to be quite robust with the interaction strength. In- 
deed, the timescale deduced from the ratio between the 
Thomas Fermi radius and the sound velocity turns out 
to be independent of gi and on the order of uJq (this 
latter result is even valid in all space dimensions and ex- 
plains why the surface modes of a dilute condensate are 
independent of the interaction strength [231 ) ■ 

The slope (parameter 7) limits ultimately the mini- 
mum achievable asymptotic velocity dispersion as illus- 
trated in Fig. [5] The value Aw(oo) obtained for t very 
small (< 3 ms) is always larger than the initial dispersion 
velocity as discussed already (see Fig. 2]). The shape of 
those curves corresponding to different values of 7 is well 




FIG. 5: Asymptotic longitudinal velocity dispersion ratio 
Au(oo)/Ai)(0) as a function of the characteristic time r used 
for the outcoupling process with an exponential time function 
f(t) = exp(— i/r): 7=1 m.s~^ (open diamonds), 7 — 0.2 
m.s~^ (open squares) and 7 = 0.01 m.s^^ (open circles). 
Same initial conditions as in Fig. [l] The dotted lines cor- 
respond to the lower bound estimated with Eq. ([SJ. 



understood: the initial decrease is due to a beneficial adi- 
abatic opening and the increase to the stretching of the 
wave function and thus to the slope. In the absence of a 
slope, only the decrease is observed. 

In order to estimate the lowest achievable asymptotic 
dispersion we proceed in the following way. For each 
value of the slope coefficient 7, we calculate with a gaus- 
sian ansatz the maximum radius size R of the wavefunc- 
tion for the lowest depth compatible with a stationary 
solution. By decreasing further the trap depth, the ab- 
sence of a stationary solution corresponds to the onset of 
the outcoupling process, and thus to the breaking of the 
adiabaticity. The corresponding asymptotic velocity dis- 
persion is then obtained by energy conservation at this 
threshold value for the Thomas Fermi radius: 
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(5) 



where Ki{R) is the proper kinetic energy per particle of 
the BEC in the Thomas Fermi regime (see Appendix) 
and E™^(R) is the corresponding interaction energy per 
particle. We have represented these limits as dotted lines 
in Fig. [5] These minima turn out to be robust against 
the specific shape of the time function f{t). They give 
a good estimate of the minimum achievable asymptotic 
velocity dispersion. For a residual slope parameter of 
7 = 0.01 m.s~^, the achievable velocity dispersion is ten 
times lower than the recoil velocity. Such a technique has 
been used in [S] to prepare a BEC for probing an optical 
lattice. 




trap depth. In ID, this relation is given by 
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FIG. 6: Output flux as a function of time for a time evolu- 
tion function (i) f{t) — exp(— f/r) (dotted line) and (ii) f{t) 
(solid line) deduced from Eq. (|6]) with a desired square flux 
Jo (dashed line). Same initial parameters as in Fig. [3] 



Kicpth(i) = M°/(0 = / 



iV(0) 



Mt)dt 



-I 2/3 



(6) 

for t > td- Figure [5] provides an example of flux engi- 
neering. Jo is taken with a square (dashed line) form in 
order to impose a nearly constant flux for a large amount 
of time. The output flux (solid line) obtained from the 
numerical simulation performed with Vdcpth(i) given by 
Eq. ([6]) gives a satisfying result compared to the desired 
flux. 

The function f{t) has also a strong influence on the 
asymptotic velocity dispersion once all atoms are out- 
coupled. As an example, for the square prescription Jq, 
the constraint on the atomic flux yields an increase by 
a factor five of the asymptotic velocity dispersion com- 
pared to the optimal two stage strategy. 



III. FLUX ENGINEERING 



IV. 2D OUTCOUPLING OF AN INTERACTING 
BOSE-EINSTEIN CONDENSATE 



Another important question from the experimental 
point of view is the possibility of controlling the out- 
put flux by a proper choice of the time function /(i). 
Let us consider the examples depicted in Fig. [B] With 
the same initial parameters as those of Fig. SI the de- 
crease of the trap depth is sufficiently small to allow out- 
coupling above a time t^ ~80 ms. With an exponential 
form f{t) = exp(— i/r) for the time evolution of the trap 
depth, the output flux Jit) = -AN/dt (where N{t) is 
the number of trapped atoms) obtained numerically has 
the shape of the dotted curve. 

In the following, we model the potential as a linear 
slope plus a harmonic potential whose angular frequency 
is kept constant. This simplification allows for an ana- 
lytic solution and does not affect the generality of our 
approach and its conclusion. We show how a quasi- 
equilibrium assumption helps to engineer the outcoupled 
atomic flux. 

The outcoupling is provided here by reducing the depth 
of the harmonic potential but not its strength. For 
t < td, the chemical potential /.t^ is smaller than the trap 
depth VdGpth(i < td) and no atoms are outcoupled. For 
t > td, we assume that the timescale over which the trap 
is decreased is sufficiently low to ensure the locking of 
the chemical potential to the instantaneous trap depth: 
M(i) = Vdepth{t > td). Assuming that the relation be- 
tween the chemical potential and the number of atoms is 
still valid in this quasi-equilibrium regime, we can infer 
the time evolution of the number of atoms N{t) and thus 
of the atomic flux J{t). 

Conversely, if we impose a desired flux function Jo{t), 
we can provide the corresponding time-evolution of the 



The ID models developed in the Sec. I and II give the 
optimal strategy for minimizing the longitudinal veloc- 
ity dispersion. In this section we consider a 2D model 
which therefore accounts also for the transverse degrees 
of freedom. Our starting point is the 2D Gross-Pitaevskii 



ih 



~dt 



2m 






dhp_ 



Uix,y,t)i^ + g2\^f^, 
(7) 



where the wave function -0 is normalized to the total 
number of particles. A first set of simulations has been 
performed with the 2D trapping potential 



Uiix,y,t) = -mjx - C/o/(i)e-2-V-„^ + 



2 2 



(8) 



where the guide is accounted for by a transverse harmonic 
potential of angular frequency ujy and /(t) = e~*''^. In 
practice, the Gaussian potential associated with the trap 
can be realized using a cylindrical red detuned off res- 
onance Gaussian beam with a large waist along the di- 
rection perpendicular to the propagation direction x. As 
expected, the results that we have obtained for the lon- 
gitudinal degrees of freedom are very close to those de- 
tailed in the previous section. Here, a sudden switch off 
is realized by setting Uq to zero while keeping the guide 
confinement unchanged. The asymptotic value of the lon- 
gitudinal momentum dispersion for a sudden switch off 
can be deduced from Eq. (jA7p of appendix A. 

Figure [7] gives examples of the variation of the trans- 
verse dispersion Ay as a function of time in the course 
and after the outcoupling for different characteristic 
times r = 1, 5 and 20 ms. The strength of the inter- 
actions, Ng2 = 267 h^/rri, has been chosen so to have 
the same Thomas Fermi radius as for the ID case with 
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FIG. 7: Transverse quadratic size as a function of time in the 
course of the outcoupling process for different outcoupling 
characteristic time r = 1 ms (dotted line), 5 ms (dashed line) 
and 20 ms (solid line). The gray straight line corresponds to 
the quadratic size of the harmonic oscillator associated with 
the transverse confinement of the guide. 



Ngi = 53/i^/^Wq m~^/^. As expected, a too fast out- 
coupling is accompanied with large transverse oscilla- 
tions. The longer the outcoupling time r, the lower the 
transverse excitation. A more quantitative analysis is 
summarized in Fig. [S] We have plotted an effective max- 
imum quantum number Anmax associated with the con- 
tamination of the transverse degrees of freedom using the 
formula 



(Ay) 



2 

max 



Ann 



1 



(9) 



where (A2/)niax corresponds to the asymptotic maximum 
amplitude of the transverse dispersion. Anmax thus mea- 
sures the contamination of the transverse excited states 
due to the outcoupling process. We observe a situation 
in which the first excited state population remains below 
10 % for T > 10 ms for our parameters. Note, however, 
that the amount of contamination for a given r increases 
slightly with the strength 32 of the interactions. 

In the limit of a sudden release of the condensate into 
the waveguide, the expansion can be accounted by a scal- 
ing law approach and the modal decomposition can be 
worked out [2^, [23] . One finds analytically and numeri- 
cally that a frequency mismatch between the initial trap 
and the guide angular frequency triggers transverse os- 
cillations. The larger the guide frequency the larger the 
transverse amplitude. However, if one performs a pro- 
gressive spilling of the atoms into the guide, these oscil- 
lations can be dramatically reduced similarly as in Fig.[7l 
If the guide frequency is larger than the trap frequency, 
a relatively low characteristic time, r, is required. In the 
opposite limit, a relatively large r is necessary to damp 
out the transverse oscillations. 

The potential Ui minimizes by construction the cou- 
pling between the transverse and longitudinal degrees of 
freedom. We have also investigated the same problem 



FIG. 8: Upper bound of the contamination of the transverse 
degree resulting from the outcoupling as a function of the 
outcoupling characteristic time r. Calculations are performed 
for a strength parameter of Ng2 = 267 fi /m and potentials 
Ui (open squares and diamonds) and Uu (black squares). An 
exponential decay f{t) = exp(— t/r) is chosen for open and 
black squares while a linear decay f{t < r) = l — t/r has been 
chosen for data represented by open diamonds. The potential 
Uii which couples strongly the longitudinal and transverse 
degrees of freedom turns out to give similar result as Ui when 
the outcoupling time is sufficiently long (r > 10 ms for our 
parameters) . 



with the potential 



Uii{x, y, t) — —mjx — Uoe 



't/r^-2ix' 



2 ^ 

_ao) 

closer to the one used in the experiments of Refs. [14l[l5|. 
For both potentials, the upper value An,„ax remains 
small (< 0.1) for r > 10 ms. This is in agreement with 
the qualitative estimate ( An'"* < 0.1) of the contribu- 
tion to An of the interactions proposed in Ref. [IJ]. In 
addition, the specific shape f{t) of the decrease of the 
trapping potential does not play a crucial role on the 
transverse state populations for a sufficiently large char- 
acteristic time r as illustrated in Fig. [5] where the results 
for a linear and an exponential form oi f{t) are compared. 

The long outcoupling times chosen in the experiments 
(Til [l5| are optimal to reduce the contamination of the 
transverse excited states and explain the possibility of 
reaching the quasi monomode regime for guided atom 
lasers. However, they are not optimal from the longitu- 
dinal dispersion point of view. 

As in Sec. |TT1 we find numerically that the optimal 
and robust strategy is still the two stage strategy which 
consists in an adiabatic opening of the trap (just before 
the time for which atoms can spill into the guide) fol- 
lowed by an abrupt switch off of the trapping potential. 
Indeed, the value obtained for Anmax turns out to be ro- 
bust against the sudden switch off of the trap. As a com- 
plementary procedure, one could envision to compensate 
the local slope by an external potential to enable a more 
efficient adiabatic opening. These features are quite gen- 



eral and expected to be still valid in 3D [28|. However, as 
in ID, it is not possible to engineer the atomic flux and 
to benefit from the lowest achievable asymptotic velocity 
dispersion. 

In conclusion, we have presented a detailed study of 
the spilling outcoupling mechanism to produce a guided 
atom laser in both position and momentum space. The 
main features and limits in velocity dispersion have been 
discussed and can be captured by simple models. The 
detail study of the control in momentum space of an in- 
teracting condensate reveals the importance of the two 
contributions to the kinetic energy: the proper kinetic 
energy of the condensate and the one associated with 
the time-dependent transformation of the potential ex- 
perienced by the condensate. Our study enables to opti- 
mize either the atomic flux or the momentum space of a 
quantum probe realized by spilling atoms from a BEC. 
In addition, it validates the entropic analysis performed 
in [iJl since non-adiabatic and interactions effects have 
been shown to play a minor role for the chosen experi- 
mental protocol. 
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FIG. 9: Time evolution of the kinetic energy of an interact- 
ing Bose-Einstein condensate (same initial conditions as in 
Fig. |4| held in an harmonic trap of initial angular frequency 
ojo = 27r X 150 Hz when the angular frequency is weakeaned as 
a function of time according to ij{t) = 4a;o(l — t/T)/5 + iJo/5: 
Numerical simulation (solid gray line), scaling law approach 
(see Eq. (|A6|) 1 (dashed line). A very good agreement is ob- 
served for various values of t = 1, 20, 200 ms. 
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Appendix A: Momentum dispersion of a non-ideal 
condensate in a time-dependent harmonic trap 

The time evolution of the shape of a condensate in the 
Thomas Fermi regime and under a time- variation of the 
trap frequencies can be accounted exactly fo r by a set of 
coupled equations involving scaling factors [2^, [30] ■ In 
this appendix, we show that this scaling approach cap- 
tures the size evolution of the condensate but fails to 
take properly into account the evolution of the momen- 
tum dispersion. Indeed, Thomas Fermi approximation 
yields a divergency of the kinetic energy at the boundary 
of the trapped condensate. 

The explicit calculation of the kinetic energy where the 
boundary is carefully treated through a suitable expan- 
sion is detailed in Ref. [3l[ for a 3D harmonically trapped 
BEC in the Thomas Fermi regime. 

We give here explicitly an approach which predicts cor- 
rectly the mean evolution of the momentum (p^) of a har- 
monically trapped BEC when the strength of the confine- 
ment is modified as a function of time. The initial angular 
frequency is denoted uj{0) = ujq and is assumed to be the 
same in all directions. The calculation is performed at 
d = 1,2,3 dimensions. 

In order to work out the time evolution of the velocity 
dispersion we have to calculate the time-dependent mean 



kinetic energy per particle along a given axis 



(pDJt) 

2m 



2m N 



d^{f,t) 



dx 



d.T 



(Al) 



Let us introduce no(^) the exact initial density profile 
before modifying the trap strength. This profile is es- 
sentially the one of an inverted parabola except in the 
vicinity of the boundaries [31| . Its initial Thomas Fermi 
radius i?° is related to the chemical potential ^° by 
i?° — (2p^^lmijj^)^l'^ . The kinetic energy per particle for 
a BEC at equilibrium in the harmonic trap is a function 
ofi?0: 



KA^W 






In 



(R'^ 



0^(2-d)/4 



f mgdN \ 



1/4 



B, 



(A2) 

where the dimensionless coefficients Ad, Bd {Ai = 0.50, 
A2 = 1.17, A3 = 2.50, Bi = 1.52, B2 = 1.45, B3 = 0.80) 
have been determined numerically following the proce- 
dure detailed in Ref. [31| for the 3D case. This method 
consists in matching the asymptotic form of the wave 
function in the interior of the cloud and in the boundary 
region and has been generalized here to ID and 2D case. 
The scaling method suggests to approximate the time- 
dependent condensate wave function by 



^d{r,t) 



1 



bd{tf°ybit) 



1/2 



-Am^^ (A3) 



where the scaling parameter b(t) obeys an Ermakov-like 
equation 



b + uj^it)b 



_UJq_ 



(A4) 
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with the initial conditions 6(0) = 1 and 6(0) = 0. 

To carry out the calculation of the kinetic energy, we 
first consider the ID situation. The kinetic energy term 
(|Af[) involves the space derivative of the wave function. 
The ansatz (jA3p has two position dependent terms orig- 
inating from the modulus pi{x,t) and phase Si{x,t) of 
the wave function. Its derivative thus reads 

dMx,t) ^ f dpx{x,t) ^ ^ dSi{x,t) \ ^,s,ix,t) 

dx \ dx ' dx J 

(A5) 
The modulus square \d'ipi{x,t)/dx\'^ involves two square 
terms that are calculated with different approximation 
schemes. The contribution to the kinetic energy of the 
square term \dpi/dx\'^ requires to include the contribu- 
tion of the wings to cancel out artificial divergences, and 
is equal to Ki[b{t)Ro]. The contribution to the kinetic 
energy of the square term \pidSi/dx\'^ is finite and can 
be carried out within the Thomas Fermi approximation 
yielding a term that scales as r. Such a procedure can 
be generalized to higher dimenisions. We finally get 



(pDit) Xrf[6(i)i?0] 



2m 



^0 V di y 



(A6) 



where Ci = 1/5, C2 = 1/6 and C3 = 1/7, and the re- 
lation between the chemical potential and the interac- 
tion strength is ^^i = SNgi/iR^^, pi = 2Ng2/iT{Rlf, 
lbNg-i/Si:{Rlf. 



Equation (|A6[) shows the two contributions to the 
kinetic energy. Depending on the characteristic time 
over which the dimensionless scaling parameter h evolves, 
the dominant contribution to the momentum dispersion 
changes: (i) in the adiabatic limit, the time derivative 
tends to zero and the first terms accounts for the scaled 
variation of the kinetic energy (case r = 100 ms in Fig. [9]); 
(ii) for a rapid change of the confinement such as in an 
abrupt switch off the second term dominates (case t — 1 
ms in Fig. |9]). For intermediate timescale, the two con- 
tributions yields oscillations of the momentum dispersion 
(case T — 2Q ms in Fig. [3]) . We find a remarkable agree- 
ment between the numerical integration of the Gross- 
Pitaevskii equation and formula (jA6|) . 

Note that the second term of the left-hand side of 
Eq. (jA6P that accounts for the transfer of the interaction 
energy into kinetic energy in the course of the expan- 
sion can be directly obtained by combining the moment 
method (via the Ehrenfest theorem) with the scaling pa- 
rameters [li-llS]. 

For an abrupt switch off {uj{t > 0) = 0), we can extract 
from Eqs. (jA4p and (|A6|) . the complete dynamics of the 
evolution of the momentum dispersion. In particular, the 
asymptotic increase value of the momentum dispersion 
reads 



(pIXc 



2C, 



rfMS 



E't 



(t = 0) 



2m 



(A7) 



^^l 
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